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1 . Introduction 


The reoently published Co spar International Referenoe Atmosphere 
CCIRA 1986 or CIRA-86) oontalns monthly tabulations of zonal mean wind 
from 0 to 120 km [Fleming et al. , 1990] derived from a tropospheric 
climatology by Oort [1983] and use of the gradient wind approximation 
with the temperature and pressure tables derived from satellite remote 
sensing data by Barnett and Corney [1985] for the middle atmosphere 
and the mass spectrometer and incoherent scatter (MSIS-83) empirical 
model [Hedin, 1983] for the thermosphere. The MSIS-83 model was the 
predecessor of the MSIS-86 model [Hedin, 1987] model which constitutes 
the thermosphere portion of the CIRA-86 density and temperature model 
[Hedin, 1988]. 

The CIRA-86 tabulations are the successor to CIRA 1972 (CIRA-72) 
which contained wind tabulations [Groves, 1972a] based largely on 
rocket and some radar measurements. However, the CIRA-86 tabulations 
were derived independently of the prior rocket data or tabulated 
models. Monthly wind models for the upper mesosphere and lower 
thermosphere based entirely on MF (medium frequency) and meteor radar 
data have also been prepared [Miyahara et al, 1991]. 

Thermospheric wind data from satellites and ground based incoherent 
scatter radar and Fabry-Perot optical interferometers have been 
oombined [Hedin et ail. , 1991] to generate an analytio empirical 
horizontal wind model (HVM90) of winds above 100 km using a limited 
set of vector spherical harmonios to describe spatial h-tiH temporal 
variations in the exosphere and at selected altitude nodes with cubic 
spline interpolation between nodes. The formulation of this wind 
model is analogous to the MSIS-86 density and temperature model and 
allows the user to obtain atmospheric parameters at an arbitrary 

MSIS-86 model was reoently extended [Hedin, 
1991] into the lower atmosphere to provide a single analytic reference 
model LMSISE— 90] of temperature and density from the ground to the 
exosphere. While the lower atmosphere portion was essentially derived 
from the CIRA-86 tabulations, the structure in the upper mesosphere 
was adjusted to best fit historical rocket data as well as maintain 
overall hydrostatio equilibrium while smoothly Joining the previous 
upper thermosphere model. 

It 16 the goal of the HWM93 model described herein to extend the 
formulation of the HWM90 wind model into the mesosphere to the 
surface so as to provide a description of the average (climatological) 
wind system throughout the atmosphere. The model is based not only on 
the CIRA-86 tabulations, but selected historical rocket data, previous 
rocket data based tabulations, meteor radar and MF radar data, and 
lower thermosphere incoherent soatter data previously used for HWM90. 
The new model thus represents a compromise between data souroes in the 
upper stratosphere, mesosphere, and lower thermosphere, while closely 
following CIRA-86 in the lower stratosphere and troposphere HWM90 
in the thermosphere. Model results and data comparisons are given with 
emphasis on the mesosphere and lower thermosphere and include also 
gradient winds calculated from MSISE-90. The HWM90 model parameters 
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were changed at 100 km to provide a smoother transition into the 
mesosphere hut otherwise remain unchanged in the thermosphere. 

2. Data Souroes 


The data used to generate this model were derived from published 
tabulations, figures, and original data bases as summarized in Table 
1. The number and diversity of data souroes is greatest in the 
mesosphere and lower thermosphere consistent with availability and the 
emphasis of this paper. 


In the mesosphere and lower thermosphere the techniques and sources 
represented Include incoherent soatter radar, MF radar, meteor radar, 
rooketsondes , rocket grenade soundings, and gradient winds. These data 
were given the most weight in deriving the model. Although MF radar 
data appear in later plots at their nominal attributed altitude, they 
were given considerably less weight in the fit above 100 km because 
their real and apparent altitudes are different near the total 
reflection height. In addition, tabulations in CIRA-72 [Groves, 1972a] 
GROVBS-69 [Groves, 1969], and CAO-83 (Central Aerological Observatory)’ 
[Koshelkov , 1983] summarizing largely older or different rocket data and 
tabulations from CIRA-86 [Fleming et al., 1990], providing global 
ooverage largely based on gradient wind calculations, were included with 
®PP^° x J“ atel y » faotor of two less weight. Gradient winds derived from 
MSISE-90 were weighted a faotor of ten less than the rocket and radar 
data and thus are only Included for comparison purposes below 85 km 
However, the CIRA-86 winds were weighted less above 85 km by a faotor of 
sixteen in order to give some precedence to the newer MSISE-90 
representation of the mesosphere/thermosphere transition region. 


In the stratosphere the data include rocketsondes, and rocket 
5oToS d «« SOU11<lin ^ s * CIRA “?2, CAO-83 , tabulations from CIRA-86 and 

winds (Table D- *or the troposphere, only the 
CIRA-86 tabulations and MSISE-90 gradient winds are used. 


While the data souroes are largely independent, 
significant overlaps which should be kept in mind, 
now be described with this aspeot in mind as far as 


there are also 
The data 6ets will 
is known. 


J 1 ^^® ©seller CIRA tabulations, the CIRA— 86 reference atmosphere 
winds in the stratosphere and mesosphere are inferred from satellite 
remote sounding pressure and temperature data and a thermospheric 
density and temperature model and make no direct use of rocket or 
radar wind data. The CIRA-86 winds were largely derived from the 
gradient wind equation. However, equatorial winds were based on the 
seoond derivative of pressure where the geostrophio formulation fails 
troposphere winds were based on a published data summary, an d very 
high latitude winds were based on a dynamical constraint (see Fleming 
et al., 1990). For longitudinal variations the CIRA-86 tabulations 
were supplemented with gradient wind calculations of wave 1 amplitudes 
and phases for the zonal and meridional wind [Fleming et al., 1988]. 
These gradient winds were favorably compared with radar derived winds 
by Manson et al. [1991]. Gradient winds (including equatorial winds 
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using the seoond derivative formulation) were also derived from the 
MSISB-90 density and temperature model and will be similar to the 
CIRA-86 winds, slnoe MSISB-90 is based heavily on the same CIRA-86 
pressure and temperature tabulations. However, there are differences, 
particularly in the upper mesosphere /lower thermosphere, and these are 
greatest near the equator. Also, winds derived from the MSISE-90 
model will have smoother variations as a function of latitude and 
month than the CIRA-86 tabulations since MSISE-90 represents a 
smoothed version of density and temperature variations. 

The CIRA-72 zonal winds (no meridional winds were published) are 
largely determined by rooket wind data and early radio-meteor results 
aT&llablo at that time. Much of the data was from reports and private 
souroes not readily available today, but there is some overlap with 
rooket data used here. Below 60 km there are separate tables for 
Amerioan and European longitudes and these were arbitrarily taken to 
apply to 90V and 20E. 

The CAO— 83 southern hemisphere reference tabulations (only zonal 
wind) are derived largely from Russian rocket data but probably 
overlap slightly the data used in CIRA-72. A later version 
[Koshelkov, 1990] incorporates more rocketsonde data but also radar 
data that is separately included here. 

The GROVES-69 tabulations are an older version of CIRA-72 for both 
zonal and meridional winds and include some of the early rocket data. 
Only the meridional wind from these tabulations is used here since 

CIRA-72 superseded the zonal winds but did not include meridional 
wind. 

Rocketsonde data from the Meteorological Rocket Network (MRN) were 
obtained from the NASA/Wallop6 database, which is similar to the World 
av 7?? ter format available from the National Climate Data Center at 
Asheville, NC. The data oover the time period from 1969 to 1991 
were separated into falling sphere data iSohmldlin, 1985], which make 
useful measurements to nearly 100 km, and paraohute/datasonde 
measurements which were limited to 75 km. Based on the available time 
period, only slight overlap is possible with CIRA-72. The MRN 
oonoentrated on taking data near looal noon, but data are available 
for all parts of the day for most stations. Por eaoh station, the 
data at two kilometer intervals were separated into two hour looal 
time groups. Monthly averages, determined by summing over all 
available years, were formed separately for the twelve looal time 
groups. These averages, which provide as equitable a looal time 
coverage as possible, were used as the rocketsonde input to the model 
and data comparisons in order to mini m ize tidal bias from this souroe. 

The rooket grenade data from 1960-1972 were partly Included in 
GROVES-69 and CIRA-72, but are not inoluded in the rocketsonde data. 

The more reoent MF and Meteor radar data Included here were not 
available for CIRA-72 and are presented only as supplementary data for 
CIRA 1986 [Manson, 1990]. 
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3. Model Formulation 


The HVM93 model is an extension of the HVM87 and HWM90 wind models 
[Hedln et al. , 1968; 1991] summarizing wind measurements in the 
thermosphere. Spatial (latitude and longitude) variations in the 
horizontal wind veotor are represented by an expansion in vector 
spherioal harmonios [Morse and Feshbaoh, 1953] with each expansion 
coefficient represented by a Fourier series in day of year for annual 
and semiannual variations. The expansion involves two orthogonal 
veotor fields, the divergence B field and the rotational C field. 

Solar activity and magnetic activity variations are not included below 
100 km. Hemispheric differences are allowed but are very limited. 

This report will treat only non-tidal variations, although tidal and 
non-tldal variations were fit in an Iterative fashion to produoe as 
self -consistent an overall model as possible. Quasi-biennial 
variations are not included. Zonally averaged meridional windB were 
not modelled below 45 km. Stationary wave 1 longitude variations were 
limited to 7 to 90 km. Only the rotational (C) field was used to 
represent longitude variations since the winds are nearly geostrophic, 
and thus nearly non-divergent , and it seemed unlikely that departures 
from a curl field could be extracted from the current limited data. 

Below 100 km the wind profiles are represented by a cubic spline , 
defined by oubio polynomials between specified nodes with first and 
seoond derivatives continuous across interior nodes. The nodes were 
Chosen to be at 100, 90, 82.5, 75, 07.5. 60. 52.5, 45, 37.5, 30. 22.5, 
15, 7.5, and 0 km providing a convenient division into intervals of 
approximately one scale height. The wind magnitude and altitude 
gradient are matched at 100 km with the thermosphere values, and in 
addition the altitude gradient is specified (fit) at 100 km. 

The harmonic expansion at each altitude node is limited to low 
order terms as summarized in Tables 2a and 2b, thus smoothing the 
model output in spaoe and time. The classification into symmetrical 
and asymmetrical is with respeot to reflection about the equator with 
symmetrical meaning the veotor spherioal harmonic term provides zonal 
winds which have the same direction across the equator while the 
meridional wind changes direction. The oolumn value 'm' refers to the 
longitudinal (or looal time) harmonio oontent (0 means no longitudinal 
variation, 1 the first harmonio, eto.). The 'n' value is the latitude 
harmonio order and is always equal to or larger than 'm' . if the n-m 
value is even, then the B field term is symmetric and the C field term 
is asymmetrlo. The higher the order 'n' the greater the latitude 
variability that oan be represented. Terms of order higher than those 
in the table were not found to be significant in fitting the present 
data set. In Table 2 a dash (-) means this term is not Included for 
this node. 

The determination of the harmonio coefficients for the various 
nodes of the wind profile is accomplished by a least squares fit to 
selected subsets of the data. The node to node variations of the 
harmonio coefficients were smoothed by refitting with the sum of the 
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squares of the differences between adjaoent node ooefflolents 
(multiplied by a oonst&nt) added to the usual sums of squares of data 
sinus model differences. The multiplicative (tension) oonstant was 
ohosen so that the sums of squares of the data residuals Increased by 
no more than one peroent. 

Root mean square deviations of the data from the model in different 
altitude regions are given in Tables 3a and 3b. The grenade and 
lnooherent scatter data tend to have the largest average departures 
because they were not smoothed or based on monthly averages. Natural 
variability is also high in the lower thermosphere where tides and 
eleotrodynamio effeots are Important. The differences between the 
model and rochet and lnooherent scatter data also reflect possible 
mismodeling of the tidal variations (model described in a separate 
report), although this is not a major factor below 100 km. 

4. Model Example6/ComparlBons/Di6cussion 

4.1 Yearly Average Latitude Variations of Zonal Wind 

The zonal and annual average zonal winds from the model are 
Illustrated in Fig . l . The mostly eastward flow in both hemispheres 
peaks at 30 m/ s at mid-latitudes in the southern hemisphere 
stratosphere (15 m/s in northern hemisphere). An equatorial zone of 
westward flow peaks near 15 m/ s in the stratosphere with small zones 
near the surface and upper mesosphere. Eastward winds in the lower 
thermosphere are not as large as in CIRA-86. 

Further examples of model winds and comparisons with data are shown 
in Figs. 2 and 3. Here data have been selected for rather broad 
altitude or latitude intervals and the model wind plotted vs latitude 
or altitude for the midpoint of the selected intervals. The example 
model plot will represent the model prediction at the exact altitude 
or latitude of individual measurements with a degree of faithfulness 
that depends on how rapidly the model varies with altitude, latitude, 
etc. The difference between the plotted points and the model line 
represents the average difference of each measurement , separated by 
souroe as specified in Table 1 , from the model (calculated exactly for 
that looation) and the error bars represent the one standard deviation 
soatter of this difference within the plot bin limits. There are up 
to three plots for each situation, separating data into three groups: 
gradient winds, one as published in CIRA-86 and one as derived form 
MSISB-90 ; rocket data, CIRA-72, and CAO-83 (primarily 
stratosphere/mesosphere ) and lnooherent soatter data (thermosphere) ; 
and meteor and MF radars (monthly averages) . 

In the stratosphere the winds are well defined by gradient winds, 
rocket data, CIRA— 72 and CAO— 83 (Figs. 2 fluid 3). The gradient winds 
are in good agreement with each other exoept at the equator , thus 
confirming the representation of temperature and density gradients in 
MSISE-90 . At low latitudes (Figs. 2a and 2b) the gradient winds may 
differ by 6-10 m/s, but neither is systematically in better agreement 
w ith the rocket data. Equatorial differences are not surprising given 


5 



tlis dependenoe of the derived zonal wind on the seoond derivative of 
pressure as a funotion of latitude, rather than the first, and the 
small magnitudes of the terms involved [Fleming and Chandra, 1989]. 
The quasi-biennial variation near the equator will also introduce 
some soatter. At high latitudes there is a systematic tendency for 
the rooket data, CIRA— 72, and CAO— 83 to be higher than the gradient 
winds in the southern hemisphere by 5 to 10 m/s and similarly lower 
than the gradient winds in the northern hemisphere. 


In the mesosphere the winds are fairly well defined by gradient 
winds, rocket data, CIRA— 72, CAO— 83, and Meteor/MF radars, but with 
increasing soatter and discrepancies toward higher altitudes at 
lower latitudes. The Meteor and MF radar values generally differ from 
the model by less than 5 m/s. At low latitudes, CIRA-86 is fairly 
close to the model but MSISE-90 is higher by 10 m/s (Fig. 3b). Rocket 
data show stronger eastward flow than the model in the lower 
mesosphere (closer to MSISE-90) and stronger westward flow than the 
model in the upper mesosphere (closer to CIRA-86). At midlatitudes, 
the rocket data, CIRA-72, and CAO-83 show stronger eastward flow by 5 
to 10 m/s than the gradient or meteor and MF winds (Fig. 2c). 


la ^e lower thermosphere there is considerable scatter between 
data points and set6 such that consistency is often poor. The spread 
between gradient winds is the order of 10-20 m/s. They follow the 
model in only a very rough way. The Meteor and MF radar data are 
generally clustered about the model within 5 m/s. At low latitudes 
rocket data, CIRA-72, and IS radar are mostly less than model while’ 
t m •K adi f at ^ wi ? d s are much above (Fig. 2d). At high latitudes, the 
to be below^odel ** above the model whil e the Meteor/MF radar tends 
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Yearly Average Latitude Variations of Meridional Wind 


Z ? n ^ L “ d 411111141 average meridional winds (Fig. l) are mostly 
northward in the northern hemisphere and southward in the southern 
hemisphere peaking near 3 m/s in the mesosphere and 6 m/s in the lower 
thermosphere with small regions of reverse flow near the mesopause at 

i«?S r +i at l tud ? S and the strat °P 4 use at high latitudes. The model uses 
°~;y t4e tw0 lowest symmetric harmonios so there is no hemispheric 
difference ezoept in direction of flow. v 

Model and data comparison plots similar to those for zonal winds 
are shown in Figs. 4 and 5, exoept there is no plot for gradient winds 
and the only rocket model is GROVBS-69. B 


There is considerable soatter among data points and data sets (Figs 
5 l' Consistency is often poor, arguing against the use of 
higher harmonios for a more complicated pattern. However, the rocket 
data (Fig. 5b) suggest that an alternating pattern of north/ south 
oells near 75 km in equatorial latitudes should probably be stronger 
than modeled. While the rooket data is fairly consistent in the lower 
mesosphere, Groves-69 is inconsistent with meteor and MF radar in the 
southern hemisphere. Meteor and MF radar data generally cluster about 
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the model within 5 m/s. There is up to 10 m/s differences among 
lnooherent soatter data in the lover thermosphere. 

4.3 Annual and Semiannual Variations of Zonal Wind 

The annual variation of the zonal wind has a winter eastward 
m aximum phase (Fig. 6) peaking at midlatitudes Just above the 
6tratopause with amplitudes of 70 m/s (southern hemisphere) and 60 m/s 
(northern). An annual variation with a summer eastward maximum phase 
exists near the mesopause with an amplitude of 10 m/s. The semiannual 
variation (Fig. 7) has an equinox eastward maximum at low latitudes 
with an amplitude of SO m/s in upper stratosphere and a westward 
ma ximum at the equinoxes in upper mesosphere with an amplitude of 15 
m/s. These features have been well documented in the literature (e.g. 
Angell and Korshover [19703; Groves, 1972b; Belmont [1985] and 
references therein) . Example latitudinal cross-sections are shown in 
Fig. 8 for four different months. 

In the stratosphere the annual and semiannual variations are rather 
well defined by gradient winds, rocket data, CIRA-72, and CAO-83. 
Comparison plots versus day of year for six broad latitude groupings 
are shown in Fig. 9. There are sometimes considerable differences 
such as at high southern latitudes where CIRA-72 and CAO-83 winds are 
well above the model and gradient winds during the September equinox 
(Fig. 9a) . 

In the mesosphere the annual and semiannual variations are fairly 
well defined by gradient wind, rocket data, CIRA-72, and CAO-83, as 
well as meteor and MF radar data. However, at low latitudes the 
gradient wind variation is qualitatively different from rocket and 
radar in the upper mesosphere. Reoently reported equatorial data 
[Vincent and Lesicar, 1991; Frltts and Isler, 1992] support a strong 
vestward wind in the spring equinox which is quite different from the 
gradient winds (Figs. 91 and 9m). The model semiannual amplitude 
reaches a deeper minimum near 65 km which is more in acoord with 
rocket data than with the amplitude suggested by gradient winds 
[Fleming and Chandra, 1989] . At northern mid-latitudes near 80 km 
(Fig. 9n), there are small but striking differences in the annu al 
variation observed by different techniques. Meteor radars have a 
weaker annual variation than either rocket or gradient winds, while 
the MF radars in this grouping (Saskatoon and Urbana) have an annual 
variation similar to the gradient winds and larger than the variation 
described by rocket data. The small average eastward flow in December 
from the meteor radars leads to an extremely weak (compared to the 
southern hemisphere winter) eastward mesospheric Jet during northern 
winter in the Mlyahara et al. [1991] model based only on radar data. 
Differences in the height of the summer reversal from westward to 
eastward flow [Manson et al., 1990] also contribute to differences in 
the annual variation in the 80 to 90 km region. For example. 

Saskatoon has a reversal height near 90 km and Atlanta and Kyoto have 
a reversal height near 80 km. In the lower mesosphere there is 
considerable separation between CIRA-72 and CAO-83 at high southern 
latitudes (Figs. 9a and 9b) where data has always been sparse. 
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In the lower thermosphere there are considerable differences between 
the gradient winds. Consistency Is often poor between data sets. The 
gradient winds are closer to eaoh other and the model during solstices 
than equinoxes (e.g. figs. 9f and 9p). The relatively large eastward 
winds of CIRA— 86 relative to radar data have been dlsoussed by Hanson 
et al. [1991]. Here the temperature and pressure gradients In the 
revised MSISE-90 appear to be an Improvement over CIRA-86/MSIS-86. At 
low latitudes, CIRA-86 Is particularly high relative to rocket, 
meteor, and HF radar during equinoxes (Figs. 91 and 9m). At high 
latitudes. Incoherent soatter shows larger variations In the 

northern hemisphere (June maximum) than meteor and HF radar or 
gradient winds (Fig. 9r). This may be due in part to limited data or 
biasing by magnetic disturbances when data yield is best. 

4.4 Annual and Semiannual Variations of Herldlonal wind 

The annual variation of the meridional wind (Fig. 6) has a January 
northward maximum peaking at 8 m/s near the equator in the mesosphere 
and reversing to a July maximum in the lower thermosphere above 105 
km. Only the lowest harmonic in latitude is used below the 
thermosphere , given the large data soatter, producing meridional winds 
of the same direction in both hemispheres. The semiannual variation 
is very small (Fig . 7 ) with an equinox poleward maximum at 

of 2 m/ s in upper mesosphere and equatorward maximum of 1 
m/s in lower mesosphere. Only the lowest harmonio in latitude is used 
providing winds of opposite direction in each hemisphere . Example 
latitudinal cross-sections are shown in Fig. 8 for four different 
months . 

Data comparison plots are shown in Fig. 10. There is considerable 
soatter among data points and data sets. Consistency is often poor 
between data sets. Overall model trends are clearly present in the 
data, but details are not particularly robust. While fluctuations in 
meridional and zonal wind measurements are similar (Table 3), the 
desired average signal is muoh smaller relative to the fluctuations 
for the meridional wind and thus incomplete ooverage in the rooket and 
Incoherent soatter data Is more noticeable. Seasonal patterns are 
often similar to Groves [1969] but of lesser magnitude (e.g. Fig. 

10a). Sometimes higher amplitudes would be consistent with limited 
radar data (Fig. 10c) but not so obviously consistent with similar 
data in the other hemisphere (Fig. 101). Latitudinal patterns are 
also similar to those of Haiyahara et al. [1991], although less 
detailed. In particular, our analysis did not find a region of winter 
to summer flow near 80 km at equatorial latitudes (imbedded in the 
more typical summer to winter flow) . The characteristics of the 
summer to winter flow as found in radar data and the relation to 
measured momentum deposition by gravity waves has been discussed by 
Hanson et al. [1991]. The reversal to a winter to summer flow in the 
lower thermosphere is clearly present in the HF/Heteor radar data, as 
it was in the Incoherent soatter data [Hedln et al., 1991], but again 
with station to station differences in reversal height as in the zonal 
wind. 
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4.5 Longitude Variations 

The stationary wave 1 (first harmonic in longitude) amplitudes are 
illustrated in Figs. 11 and 12. The zonal variations peak in the upper 
stratosphere at northern winter mid-latitudes with an amplitude of 30 
m/s (zonal) and again at the pole with an amplitude of 35 m/s (zonal and 
meridional) . Meridional vind variations peak at the pole with the same 
amplitude. Southern hemisphere amplitudes are less than 10 m/s. Annual 
and semiannual variation in wave amplitude are included in the model. 
Longitude variations related to pi erne t ary waves are not carried above 90 
km for lack of defining data (but longitude variations are present in 
thermosphere above 130 km because of physical processes tied to the 
magnetic field geometry) or below 15 km (limited by the CIRA 
tabulations) . 

The longitude variations are derived almost entirely from gradient 
wind data. Comparisons of the gradient winds from remote satellite 
data [Fleming et al., 1988] and from MSISE-90 with the HWM93 model are 
shown in Figs. 13 and 14. The two gradient wind estimates are 
generally similar with eaoh other except near the poles where there 
was difficulty in performing the numerical differentiation from the 
tabulated satellite data and near the equator where, as for zonal 
averages, the wind depends on a second derivative of pressure. The 
month to month variability is remarkably well represented by a sum of 
annual and semiannual variations in the longitudinal harmonic 
coefficients. Other data are nominally consistent, but sire 
insufficient in longitude ooverage to define the variation. In 
particular, the CIRA-72 tables, separated by Amerloan and European 
longitudes, are reasonably consistent with the satellite data. An 
example comparing gradient winds with other data is shown in Fig. 15. 
Systematic differences are notioed between gradient and radar winds as 
seen also in the zonal average plots (Fig. 9n). 

Seoond harmonic (wave 2) variations are present [Barnett and 
Labltzke, 1990; Fleming et al. , 1988] in the stratosphere and 
mesosphere but, unlike the wave 1 variations, their month to month 
variability is more random and not as usefully represented in terms of 
a mean and annual and semiannual variations and is thus not included in 
this model. 

5 . Summary 

Reference winds from CIRA-86 combined with rocket soundings, 
incoherent scatter, MF radar, meteor radar, and older reference 
tabulations have been used to extend the HWM90 spherical harmonic wind 
model into the lower atmosphere providing a unified description of 
zonal and meridional prevailing winds from the surf aoe to the 
exosphere. While month to month details oannot be completely 
represented, mesospheric data are fit with an overall rms error of 
approximately 15 m/s and considerably better in the stratosphere. 
Comparison with rocket and radar data indicates that the model 
represents current knowledge of climatology reasonably well. 
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Tab I# 2a 


Maximum B Field Spherical Harmonic Order (n) 


Term 


B field parameter (node altitude) 

100 90 82 75 67 60 52 45 37 36 22 15 7 


grad 
0 100 


Symmetrical 
Time Indep. 0 

Semiannual 6 

Asymmetrical 
Annua I 0 


Table 2b. 


Term 


Symmetrical 
Time Indep. 0 

Annual 0 

Semiannual 0 

Longitude 1 

Longitude- 1 

annual 

Longitude- 1 

semiannual 

Asymmetrical 
Time Indep. 6 

Annua I 6 

Sem I annuo I 0 

Longitude 1 

Longitude- 1 

annua I 

Longitude- 1 

semiannual 


4 4 4 4 4 

- 2 2 2 2 


3 1111 


5 5 5 5 5 

-1111 
3 5 5 5 5 

--666 
--666 

--666 


---22 
4 6 6 6 6 

- 4 4 4 4 

--555 
--555 

--555 


2 2 2 

2 2 - 


5 5 5 5 5 

11111 

5 5 5 5 5 

6 6 6 6 6 

6 6 6 6 6 

6 6 6 6 6 


2 2 2 2 2 

6 6 6 6 6 

4 4 4 4 4 

5 5 5 5 5 

5 5 5 5 5 

5 5 5 5 5 


2 


1 


5 5 5 5 5 

1 1 - - 

3 3 - - 3 

6 6 - - 

6 6 - - 

6 6 -- 


2 2 2 2 

6 6 6 6 2 

4 - - - 

5 5 - - 

5 5 - - 

5 5-- 


Maximum C Field Spherical Harmonic Order (n) 


C field parameter (node attitude) 

grad 

m 100 90 82 75 67 60 52 45 37 30 22 15 7 0 106 
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TABLE 3a. RMS Zoned Wind Differences from HWM92 


Data Set 

15 

1 o 
1 <0 
i i 

1 o 
1 <0 
1 
1 

-90 km 

90- 

120 km 


rms 

pts 

rms 

pts 

rms 

pts 

MSIS Gradient 

8 

1428 

12 

1020 

15 

1020 

CIRA-86 

8 

2521 

13 

1973 

21 

1825 

Datasonde 

7 

4949 

13 

1694 



Falling Sphere 

10 

1550 

13 

1336 

32 

135 

Grenade 

16 

771 

18 

973 

55 

146 

CIRA-72 

16 

1655 

18 

1225 

22 

744 

CAO-83 

10 

768 

12 

480 



MF radar 



6 

1029 

6 

896 

Meteor radar 



11 

776 

6 

866 

IS radar 
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TABLE 3b. RMS 

Meridional Wind Differences from HWl! 

Data Set 

16- 

A j 

O 1 
<0 1 

60- 

I 1 

1 o 
i a> 

90-120 km 


rms 

pts 

rms 

pts 

rms 

pts 


Datasonde 
Falling Sphere 
Grenade 
GROVES-69 
MP radar 
Meteor radar 
IS radar 


4 4949 

4 1551 

10 772 

6 204 


9 1684 

9 1368 

16 973 

8 1260 
5 1028 

9 394 


27 138 

50 147 

12 972 

5 896 

6 487 

34 2064 


r0 ? t mean s S uare difference between data and model, pts is 
soatter° f sample P olntB> ** is medium frequency, and IS is incoherent 
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SDL AVERAGE ZONAL WIND 




- 6 to 

1 o- 


LAT I TUDE 

SDL AVERAGE MERIDIONAL WIND 



LATITUDE 


Pig. 1. Contour plots In altitude versus latlt _ 

diurnal, and longitudinal (SDL) average zonal wind (upper panel) '‘and 
meridional wind (lower panel). Positive is eastward and northward 
(negative dashed) . H and L indicate relative highs and lows 
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SDL Average 
25 - 35 ka 


lira 


SDL Average 


LATITUDE 


LATITUDE 



LATITUDE 



LATITUDE 


MBS 


LATITUDE 


LATITUDE 


£*2- 2a. SDL average zonal wind versus latitude for 25 to 35 and 35 

Piot 6 fiS?hn? 5 e H 5 I? 2 3 + Wl ? d 4 . (s0lid lin ®) shown for mid-range conditions 
Plot symbois lndloate data souroe as given in Table 1. Top row of 

winds , middle row rocket data, and Bottom row 
MF /Meteor radar data. See text for further description. 












z ° n ® >1 Wind versus latitude for 45 to 55 and 55 
£°® 5kia * The HWM93 wind (solid line) shown for mid-range conditions. 
Plot symbols lndioate data souroe as given In Table 1. Top row of 

2 r 2 d i ent winds, middle row rooket data, and. bottom row 
MF /Meteor radar data. See text for further description. 
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fi g inR d v« S SJL a 2S£5l e versus latitude for 85 to 95 and 95 

p?ni°BTri£Ai F 1 ? SJ 119 ? w i n £ (solid line) shown for mid-range conditions. 
Plot symbols Indicate data souroe as given In Table 1 . Top row of 

ISSJ JL°22 t# w5 /I r ? dlent winds . middle row rocket and IS data. SSd 
bottom row MF/leteor radar data. See text for further description. 
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MERIDIONAL WIND (a/sec) MERIDIONAL WIND (m/sec) 


Pig. 6a. SDL average meridional wind versus altitude for southern 
high and middle latitudes. The HWM93 wind (solid line) shown for 
mid-range conditions. Plot symbols indicate data source as given ii 
Table 1. Left column of plots contains rochet data and right oolumi 
MF /Meteor radar data. 6 






SDL 

Ave rage 
30S - 0 


-20 - 1 0 0 10 20 

MERIDIONAL WIND (a/sec) 



-20 -10 0 10 20 
MERIDIONAL WIND (m/sec) 



-n-tJt.— 


SDL 

Ave rage 


-20 - 1 0 0 10 20 
MERIDIONAL WIND (m/sec) 



UJ 

§ 60 



-20 -10 0 10 20 
MERIDIONAL WIND (m/sec) 


Til;* 6 ?' S ££ meridional wind versus altitude for equatorial 

latitudes. The HWM93 wind (solid line) shown for mid-range conditions. 
JJoJ J B 2 ab 2 1 ? iadioat ® data source as given in Table 1. Left column of 
plots contains rochet and IS data data and right column MF/Meteor 
r&Cl&r Q8lt&* 
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MERIDIONAL WIND (m/sec) MERIDIONAL WIND (m/sec) 

Fig. 5c. SDL average meridional wind versus altitude for northern 
high and middle latitudes. The HWM93 wind (solid line) shown for 
mia-range conditions. Plot symbols indicate data source as given in 
Table 1. Left column of plots contains rochet and IS data and right 
column MF /Meteor radar data. 
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P L av ®J* a 2 e zonal wind versus day of year for 35 to 45 and 55 
to 65 km at southern high latitudes. The HWM93 wind (solid line) shown 
for mid range conditions . Plot symbols indioate data source as given 
2* To ? plots contain gradient winds, middle row 8 

rochet data, sad bottom row MF/Heteor radar data. 
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Pig. 9b. DL average zonal vlnd versus day of year for 65 to 75 and 75 
to 85 km at southern high latitudes. The HWM93 wind (solid line) shown 
for mid-range conditions. Plot symbols indicate data souroe as given 
in Table 1. Top row of plots oontain gradient winds, middle row 
rocket data, ana. bottom row HF/Heteor radar data. 


37 
























Fig. 9d. DL average zonal wind versus day of year for 35 to 45 and 55 
to 65 km at southern mid latitudes. The HVM93 vlnd (solid line) shown 
for mid-range conditions. Plot symbols Indicate data source as given 
In Table 1. Top row of plots oontaln gradient winds, middle row 
rocket data, ana bottom row MF /Meteor radar data. 
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102.00 

DAY OF YEAS 


• pL average zonal wind versus day of year for 65 to 75 and 75 
to 85 km at southern mid latitudes. The HWM93 wind (solid line) shown 
7 0r J n 2-? -r ? J1 2 e conditions. Plot symbols indioate data souroe as given 
in Table 1. Top row of plots contain gradient winds, middle row 
rooket data, and bottom row MF /Meteor radar data. 
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DAT OF TEAR DAT OF TEA* 

Fig. 9f . DL average zonal wind versus day of year for 85 to 95 95 

£° 1Q 5 km at southern mid latitudes. The HWM93 wind (solid line) shown 
i or «.®r, ran 2® conditions . Plot symbols indicate data source as given 
in Table 1. Top row of plots oontain gradient winds, middle row 
rocket data, ana bottom row MF/Meteor radar data. 
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* P* av ®£ a 2 e zonal wind versus day of 

latitudes . The HWM93 wma ^BOiia ±i.nej snown 
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OAT OF YEAR DAT OF TEAR 


Fig. 9b. . DL average zonal wind versus day of year for 65 to 75 and 75 
to 85 km at southern low latitudes. The HWM93 wind (solid line) shown 
for mid-range conditions. Plot symbols indicate data source as given 
in Table 1. Top row of plots contain gradient winds, middle row 
rooket data, ana bottom row MF/Meteor radar data. 
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Fig. 9j. DL average zonal wind versus day of year for 35 to 45 and 55 
to 65 rm at northern low latitudes. The HwM93 wind (solid line) shown 
for mid-range conditions. Plot symbols indicate data source as given 
in Table 1. Top row of plots oontaln gradient winds, middle row 
rochet data, ana bottom row MF/Meteor radar data. 
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pi#- 9 £- pi* average zonal wind versus day of year for 65 to 75 and 75 
to 85 km at northern low latitudes. The HVM93 wind (solid line) shown 5 

fS r *nIi?: r ? nge H. oon<LLtio 5 s *, ? lot B 7 mb °ls indicate data souroe Is gi?eS 
I' To ? P lots contain gradient winds, middle row 8 

rochet data, and bottom row MF/Meteor radar data. 
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Cl .SC IC2.CC 17C.7S CM. 00 

DAY OF YEAR 


Fig. 91. DL average zonal wind versus day of year for 85 to 95 and 95 
to 105 km at northern low latitudes. The HWM93 wind (solid line) shown 
for mid-range conditions. Plot symbols Indicate data source as given 
In Table 1. Top row of plots contain gradient winds, middle row 
rocket and IS data, and bottom row MF/Meteor radar data. 
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P L av ®J a 8® z ??al wind versus day of year for 35 to 45 and 55 
to 65 km at northern mid latitudes. The HVM93 wind (solid line) shown 
for aid-range oonditions. Plot symbols Indicate data souroe as given 
I * To ? £°w of plots contain gradient winds, middle row 8 
rocket data, and bottom row MF/Meteor radar data. 
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Pig. 9n. DL average zonal wind versus day of year for 65 to 75 and 75 
to 85 km at northern mid latitudes. The HwM93 wind (solid line) shown 
for mid-range conditions. Plot symbols indicate data source as given 
in Table 1. Top row of plots contain gradient winds, middle row 
rocket data, ana bottom row MF/Meteor radar data. 
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®°* DJf average zonal wind versus day of year for 85 to 95 and 95 
3® northern mid latitudes. The HVM93 wind (solid line 1 ) shown 

r q , l,, i? _ran g® ooaditions. Plot symbols indicate data souroeas divan 
2‘ To ? l ov of P lot B oontain gradient winds middll ?o W g 
rooket data, anfi bottom row MP/Meteor radar data row 
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BAT OF YEAR DAY OF YEAR 



Fig. 9p. DL average zonal wind versus day of year for 35 to 45 and 55 
to 65 km at northern high latitudes. The HWM93 wind (solid line) shown 
for mid-range conditions. Plot symbols indicate data source as given 
in Table 1. Top row of plots contain gradient winds, middle row 
rocket data, ana bottom row MF/Meteor radar data. 
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9 3- P L average zonal wind versus day of year for 65 to 75 and 75 
Jpes at northern high latitudes. The HWM93 wind (solid line) shown 
fS^?: r ^ oonditio $ e -, symbols indicate data source as given 

iJLwS 1 ? 1* To ? £°w of plots contain gradient winds, middle row 8 
rocket data, and bottom row MF/Meteor radar data 
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Fig. 9r. DL average zonal wind versus day of year for 85 to 95 and 95 
to 105 km at northern high latitudes. The HWM93 wind (solid line) 
shown for mid-range conditions. Plot symbols Indicate data source as 
given In Table 1. Top row of plots contain gradient winds, middle row 
rocket and IS data, and bottom row MF/Meteor radar data. 
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Fig. 101. DL average meridional wind versus day of year for 85 to 95 and 95 to 105 km at 
northern high latitudes. The HWM93 wind (solid line) shown for mid-range conditions. Plot 
symbols indicate data source as given in Table 1. Top row of plots contains rocket and IS 
data and bottom row MF/Meteor radar data. 
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Pig. 11. Contour plots in altitude versus latitude of the seasonal and diurnal (SD) 
average wave 1 amplitude and phase of the zonal and meridional wind. 
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LATITUDE LATITUOE 

Pig. 12a. Contour plots in altitude versus latitude of the wave 1 amplitude and phase for 
the zonal and meridional wind for January (midmonth). <unpxn;uae ana pnase for 
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Pig. 14a. Amplitude and phase of wave 1 zonal and meridional wind versus day of year for 
55 to 65 km at southern high latitudes. The HWM93 model (solid line) is shown for 
mid-range conditions and compared to gradient wind calculations. 
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Fig. 15. Zonal wind versus longitude for model and data at northern 
mid-latitudes and selected altitudes. 
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